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Abstract. Neutrinos are a very promising mes- 
senger at tens of EeV and above. They can be 
produced by several channels, namely as by products 
of hadronic interactions at the sources, as the main 
products of the decay of super massive particles and, 
in a guaranteed way, as the result of the propagation 
of UHECR through the bath of microwave relic 
photons. A new era of very large exposure space 
observatories, of which the JEM-EUSO mission is a 
prime example, is on the horizon and, with it, it is 
even larger the possibility of astrophysical neutrino 
detection at the highest energies. In the present work 
we use a combination of the PYTHIA interaction 
code with the CONEX shower simulation package in 
order to produce fast one-dimensional simulations 
of neutrino initiated showers in air. We make a 
detail study of the structure of the corresponding 
longitudinal profiles, but focus our physical analysis 
mainly on the development of showers at mid and 
high altitudes, where they can be an interesting target 
for space fluorescence observatories. 
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I. Introduction 

The neutrino flux carries very important astrophysical 
information. A high energy neutrino flux is expected as 
a by-product of the interactions of cosmic ray hadrons 
at the sources [1]. They can also be produced during the 
propagation of cosmic rays in the intergalactic medium 
111 and as the main product of the decay of superheavy 
relic particles |3|, |4|. 

JEM-EUSO |5| with its 10^^ tn of atmospheric target 
volume has the real possibility of observing high energy 



neutrinos and make important contributions to the under- 
standing of UHECR production and propagation Q. 
Source distributions rapidly evolving with redshift would 
be particularly favorable by increasing the cosmogenic 
neutrino flux at highest energies [2|. A thorough under- 
standing of neutrino deep inelastic scattering, as well 
as the evolution of longitudinal profiles of atmospheric 
neutrino showers, are extremely important in order to 
take advantage of the full potential of the experiment. 
Conversely, besides the obvious astrophysical value, the 
properties of just a few observed showers can also 
give valuable information on the physics governing high 
energy neutrino-nucleon interactions. The objective of 
this work is to present the first part of an ongoing effort 
in that direction. 

II. Neutrino nucleon interaction 

High energy neutrinos that propagate in the Earth 
atmosphere can interact with protons and neutrons of 
the air molecules. There are two possible channels for 
this interaction, charge and neutral current. The major 
uncertainty on the differential cross section at the en- 
ergies considered comes from the unknown behavior of 
the parton distribution functions (PDFs) at very small 
values of the parton momentum fraction x. 

The simulation of neutrino nucleon interaction is 
performed by using the PYTHIA code fS]. The parton 
distribution library LHAPDF |9| is finked with PYTHIA 
in order to be able to use different extrapolations of 
the PDFs. In this work the CTEQ6 flOl and GJR08 
1 11 1 PDF sets are considered. Fig. [T] shows the energy 
fraction carried by the leading particle as a function 
of the incident neutrino energy for the charge current 
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Fig. 2. Average energy fraction for tlie most important particles 
recognized by CONEX produced in a cliarge current interaction of an 
electron neutrino with a proton for thi'ee different neutrino energies. 

interaction of an electron neutrino with a proton for 
both sets of PDFs considered. In both cases the energy 
fraction increases steadily with the incoming neutrino 
energy. The difference between both PDFs increases up 
to a maximum of a few percent at the highest energies. 

Besides the leading particle, different types of sec- 
ondaries are generated as a result of the interaction. In 
particular, in this work we are interested in the ones 
recognized by CONEX lfT2l code which is used to 
simulate the neutrino showers. Fig. |2] shows the energy 
fraction taken by the most relevant particles recognized 
by CONEX for three different electron neutrino energies, 
obtained as a result of the charge current interaction 
with a proton. It can be seen that the smaller the 
neutrino energy the larger the energy fraction taken by 
the secondary particles. 

III. Neutrino showers 

The particles produced in a neutrino-nucleon inter- 
action are injected in CONEX with QGSJET-II 113] 
producing extensive air showers. Because the mean free 
path of neutrinos propagating in the atmosphere is very 
small, they can interact very deeply, after traversing a 
large amount of matter Fig. [5] shows the energy deposit 



as a function of X — Xq, where Xq corresponds to the 
injection point, for horizontal electron neutrino showers 
at sea level of = 10^" eV, injected at Xq = 36500 g 
cm^^ (maximum slant depth for a completely horizontal 
shower). Note that the interaction point is situated on the 
vertical axis of the JEM-EUSO field of view (FOV). Fig. 
|3] shows very broad profiles which may present several 
peaks and large fluctuations. This behavior is due to 
the Landau Pomeranchuk Migdal (LPM) effect, which 
is very important inside dense regions of the atmosphere 
and for electromagnetic particles, electrons in this case, 
which take about 80% of the parent neutrino energy. 

An orbital detector like JEM-EUSO can also detect 
horizontal showers that do not hit the ground. In partic- 
ular, horizontal neutrinos can interact at higher altitudes 
producing a shower observable by the detector. Fig. 
m shows the mean value and one sigma confidence 
level regions for the longitudinal profiles for horizontal 
electron neutrino air showers of E^, = 10^° eV, injected 
at different altitudes in points contained on the vertical 
axis of the cone of the FOV. It can be seen that as 
the altitude increases the fluctuations are reduced and, 
on average, the profiles become thinner. This is due 
to the fact that the LPM effect become progressively 
less important with increasing altitude because of the 
decrease in atmospheric density. Note that just showers 
up to 20 km of altitude are considered because, at higher 
altitudes, the grammage of the FOV is not enough to 
contain the whole profiles. 

As already mentioned, at smaller altitude the longitu- 
dinal profiles present a complicated structure (see Fig. 
|3]l. In particular, the showers present multiple peaks. 
Fig. |5] shows the distribution functions of the position 
of each maximum. X^^^. corresponds to the position 
of the first maximum counted from the start of the 
shower, X^^^^^ is the second one and so on. Note that the 
distribution function of X^^^ is bi-valued and its first 
peak is located at ~ 800 g cm^^, while the second one 
is at ~ 1500 g cm^^. The first peak corresponds to the 
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Fig. 4. Mean value and one sigma regions of the longitudinal profile 
corresponding to horizontal electron neutrino showers of Ei, = 10^" 
eV for different altitudes. The interaction point is contained on the 
vertical axis of the JEM-EUSO telescope. 
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Fig. 5. Distribution of the position of the maxima for electron neutrino 
showers of Ei, = 10^" eV at sea level. 

development of the hadronic component of the electron 
neutrino cascade, whereas the second one reflects the 
electromagnetic portion of the shower. 

This complicated structure simplifies as the altitude 
increases. Fig. |6] shows the probability of finding a 
profile with a given number of peaks, Nxi ■ As the 
altitude increases the probability of finding a shower 
with more than one peak goes to zero. 

The first portion of the cascades is, in general, dom- 
inated by the hadronic component. Therefore, the way 
in which the primary neutrino energy is distributed can 
be assessed by, 

F,,, = /o"^ '^^ ^ ^ (1) 

where dE/dX is the energy deposition, Xc — 1200 
gcm^^ is a characteristic depth that roughly sepa- 
rates the hadronic -dominated from the electromagnetic- 
dominated portions of the shower, and Xum is the 
maximum atmospheric depth reached by the shower. F^n 
is calculated only for those showers that present a first 
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Fig. 6. Probability of an electron neutrino shower to have A'^j,^; 
maxima for E^, = 10^*^ eV for horizontal showers at different 
altitudes. 

maximum at a depth smaller than Xc. Fig. |7] shows the 
Fen distributions obtained for CTEQ6, GJR08 and for a 
modification of the CTEQ6 results in which the leading 
particle takes 70% of the neutrino energy, corresponding 
to horizontal 10^" eV electron neutrinos injected at sea 
level at the axis of the FOV. It can be seen that F^n is 
correlated with the energy taken by the leading particle. 
The smallest value of the average of F^n corresponds to 
CTEQ6 because the leading particle takes about 82% of 
the neutrino showers while, in the case of GJR08, it takes 
~ 79%. The differences between the results obtained 
for these models are very small. In the case in which 
the energy taken by the leading particle is artificially 
reduced to 70% (hashed histogram), the mean value of 
Fen increases. Note that, in the latter case, the energy 
extracted from the leading particle in order to reduce its 
average energy, is redistributed among the other daughter 
particles in such a way that the ratio between the energy 
taken by a given particle and the total energy taken by all 
secondaries, excluding the leading particle, is constant. 

The sensitivity of F^n to the hadron component of 
the showers depends on the altitude. Showers injected 
at sea level in regions of high density are dominated by 
the LPM effect allowing a clearer separation between the 
hadronic and electromagnetic portions of the cascades. 

The parameter F^n can be very useful to understand 
neutrino interactions with atmospheric nuclei and, in 
particular, to estimate the energy fraction taken by the 
leading particle. Any practical application, however, will 
depend on the actual event rate. 

IV. Proton and neutrino events 

The interaction length for protons is \pr{lQ^^eN) ~ 
36 g cm-2 ^jjjj fQj. neutrinos is A,,(1020eV) - 3.2 x lO"^ 
g cm^^. The survival probability of an horizontal proton 
that reaches the Earth surface at the vertical axis of 
the FOV is ^ exp(— 1000), whereas the corresponding 
probability for a neutrino is ~ exp(— 0.001). Therefore, 
despite the fact that horizontal neutrino and proton 
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Fig. 7. Distribution function of parameter F^n (see Eq. Q) for 
CTEQ6, GJR08 models and also for an artificial modification of the 
CTEQ6 prediction such that the leading particle takes 70% of the 
neutrino energy. Horizontal showers of 10^" eV injected at sea level 
at the axis of the FOV are considered. 



showers have very different observational characteristics, 
it is very unlikely to observe a proton interacting so deep 
in the atmosphere. 

Nevertheless, for a given proton and a neutrino fluxes, 
there exists a particular slant depth for which the proton 
and neutrino events have the same rate. For any par- 
ticle type, the probability of interacting in the interval 
\X, X + AX] is given by. 



Prnt{E;X, AX) 



where A(i?) 



exp(-X/A(£;)) X 

[1 - exp(-AX/A(S))] , (2) 



is the interaction length at a 
given energy. Therefore, solving the equation 
c^priE) PZiE;Xo,AX) = ME) PUE;Xo,AX), 
where (j>pr{E) and 4)^{E) are the proton and neutrino 
fluxes, the slant depth at which protons and neutrinos 
can be detected with the same rate is. 
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exp(-AX/A,(S)) ; 

The function X{E) is obtained by using: AX ~ 
Xpr{E), the proton-air cross section of Sibyll 2.1 
lfT4l . the neutrino cross section a'^^{E) = 6.04 x 
10-^^{E/GeW)°-^^^ cm2 [15|, the Waxman-Bachall up- 
per limit for the neutrino flux | i6| and a power law fit 
of the Auger spectrum ifTTl . Therefore, for E = 10^° 
eV the event rates for protons and neutrinos are of the 
same order of magnitude for X = 142 gcm^^. The 
latter means that, under the assumptions of the present 
calculation, protons can act as a background for neutrino 
identification in inclined events. 

Fig. [8] shows the average profile and the 68% CL 
for proton and electron neutrino induce air shower 
longitudinal profiles of zenith angle 6 = 45° and primary 
energy 10^" eV injected at 142 gcm^^. Although, the 
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Proton and electron neutrino air showers of 10^" eV, 6 = 45° 
at 142 gcm~2. 



event rate of this kind of showers are similar for both 
protons and neutrinos, the profiles, however, are quite 
different making a good discrimination possible. More 
specifically, it can be seen that the neutrino showers 
develop deeper in the atmosphere and present larger 
fluctuations, due to the LPM effect. 

V. Conclusions 

Neutrino detection is of great importance for the 
understanding of several astrophysical process and, in 
particular, the origin and propagation of the highest en- 
ergy cosmic rays. JEM-EUSO is technically capable of 
observing neutrino initiated cascades and to discriminate 
these cascades from hadronic ones. Furthermore, the 
high degree of structuring of the longitudinal neutrino 
shower profiles opens a spectrum of experimental oppor- 
tunities which JEM-EUSO will be able to explore, with 
larger or lesser success, depending on the exact shape 
of the neutrino energy spectrum at the highest energies. 
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